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ABSTRACT 

Near-infrared medium-resolution spectra of seven bright brown dwarfs are presented. The 
spectra were obtained with the Infrared Camera (IRC) on board the infrared astronomical satellite 
AKARI, covering 2.5-5.0 fj,m with a spectral resolution of approximately 120. The spectral types 
of the objects range from L5 to T8, and enable us to study the spectral evolution of brown dwarfs. 
The observed spectra are in general consistent with the predictions from the previous observations 
and photospheric models; spectra of L-type dwarfs are characterized by continuum opacity from 
dust clouds in the photosphere, while very strong molecular absorption bands dominate the 
spectra in the T-type dwarfs. We find that the CO fundamental band around 4.6 /im is clearly 
seen even in the T8 dwarf 2MASS J041519— 0935, confirming the presence of non-equilibrium 
chemical state in the atmosphere. We also identify the CO2 fundamental stretching-mode band 
at 4.2 /im for the first time in the spectra of late-L and T-type brown dwarfs. 

As a preliminary step towards interpretation of the data obtained by AKARI, we analyze 
the observed spectra by comparing with the predicted ones based on the Unified Cloudy Model 
(UCM). Although overall spectral energy distributions (SEDs) can be reasonably fitted with the 
UCM, observed CO and CO2 bands in late-L and T-dwarfs are unexpectedly stronger than the 
model predictions assuming local thermodynamical equilibrium (LTE). We examine the vertical 
mixing model and find that this model explains the CO band at least partly in the T-dwarfs 
2MASS J041519-0935 and 2MASS J055919-1404. The CO fundamental band also shows excess 
absorption against the predicted one in the L9 dwarf SDSS J083008+4828. Since CO is already 
highly abundant in the upper photospheres of late-L dwarfs, the extra CO by vertical mixing has 
little effect on the CO band strengths, and the vertical mixing model cannot be applied to this 
L-dwarf. A more serious problem is that the significant enhancement of the CO2 4.2 /im band 
in both the late-L and T dwarfs cannot be explained at all by the vertical mixing model. The 
enhancement of the CO2 band remains puzzling. 

Subject headings: molecular processes — stars: atmospheres — stars: late-type — stars:low-mass, 
brown dwarfs — 
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1. Introduction 

Brown dwarfs are denned as the objects that 
were born as isolated objects but were not mas- 
sive enough to ignit e hydrogen nuclear burning 



( Burrows et al. 2001). They are particularly of 



interest for their extremely low-temperature at- 
mospheres in relation to giant planets. Observa- 
tional studies of brown dwarfs have been exten- 
sively carried out since the first definitive identifi- 
cation of a genuine brown dwarf Gl 229B (of the 
class n ow known as T-dwarf) by iNakaiima et al. 
(1 19951) . The firs t L-type brown dwarf GD 1 65B 



was reported by iBecklin fc Zuckermanl ()l988l ) al- 
though it was debated for some time whether it 
was a star or a brown dwarf. This very red object 
is important in that the presence of dust layer in 
the photospher e was recognized f or the first time 



by this object ( Tsuii et al. 19961 ). Spectroscopic 



observations in the infrared regime are the most 
essential tools to obtain physical and chemical in- 
formation of brown dwarfs, through absorption 
bands of various molecules in the atmospheres. 

Carbon monoxide (CO) plays a decisive role 
in determining the major characteristics of cool 
stars because o f its very large dissociation energy 
(jRussell I Il934j ). In brown dwarfs, the CO first 



overtone band at 2.3 /i m is observed until e arly- 
T type sources ( s ee e.g. , Geballe et al. ( 20021 ) and 
Burgasser et"afl (|2006bl) ). The role of CO, how- 
ever, changes in very cool dwarfs: carbon resides 
mostly in CH4 rather than in CO in a very cool 
(e.g. T w 1000K) and high de nsity (e.g. logP g w 
6.0) environment ( Tsuiil 19641 ) . and CO no longer 
plays any critical role under such a circumstance. 
Although no celestial object outside our Solar sys- 
tem having such a physical condition had been 
known, it was found at last that such a case is ac- 
tually realized in the b rown dwarf Gl 229B show - 



ing strong CH4 bands ( Oppenheimer et al. 19951 ) 



However, nature is not so simple as predicted 
by a simple theory of thermochemistry. In fact, it 
was not long before an unexpected detection of CO 
fundamental band at 4.6 in Gl 229B, classified 
as T6 , was reported bv |Noll. Geballe. k, Marlev 
(|l997t) and lOppenheimer et all (|l998l) . These pi- 
oneering observations highlighted roles of non- 
equilibrium chemical processes. The idea of ver- 



tical mixing_was_su£gested_by iGriffith fc Yelle 



extreme stability of CO, the chemical time scale 
is much longer than the mixing timescale and CO 
is dredged up by vertical mixing from the inner 
part of the atmosphere where CO is still abun- 
dant to the surface layers where CH4 dominates 
instead of CO. Recently two more late-T dwarfs 



were f ound to show the CO band bv I Geballe et al 
(2009), confirming that the presence of CO is a 



general characteristic of late-T dwarfs. 

It is to be remembered that the observations 
so far are limited to late-T dwarfs. Observations 
of the CO fundamental band in early-T and L- 
type dwarfs are needed for a unified understanding 
of the physical and chemical processes related to 
CO in the atmospheres of brown dwarfs. Spectro- 
scopic observations of the CH4 and of the CO tran- 
sitions in L and M-band are always difficult from 
the ground. Severe atmospheric absorption and 
limited wavelength coverage make precise analy- 
sis difficult. Thus, there is a strong motivation 
to carry out near-infrared spectroscopy of brown 
dwarfs from space, especially in the region of 2.5- 
5.0 /im which remains the least explored so far. 

The infrared astron omical satellite AKARI 
( Murakami et al.l 120071 ) was launched in 2006 



(ll999f ) and lSaumon et all (|2000h . Because of the 



February. Scientific observations under cryogenic 
condition were carried out from 2006 May to 2007 
August, and ended with the liquid Helium boil- 
off. The satellite is equipped with a 68.5 cm cooled 
telescope and two scientific instruments covering 
the wavelength range of 1. 8-180 /xm. One of t hem, 
the Infrared Camera (IRC: IOiiaka et al l l2007l) car- 
ried out imaging and spectroscopic observations 
in the near- to mid-infrared wavelength regions. 
The IRC provided a unique opportunity to take 
spectra of brown dwarfs in the range between 2.5 
and 5.0 /im without interference of telluric ab- 
sorption from the atmosphere. In the framework 
of the AKARI Mission Programme (MP; coordi- 
nated observations by the project team members), 
we have carried out a series of near-infrared spec- 
troscopic observations of ultracool dwarfs. This 
programme, NIRLT (PI: I.Yamamura), aims to 
obtain a set of legacy data for studies of the evo- 
lution of physical and chemical structure of L- 
and T-dwarfs. In this paper we present the initial 
results from this programme. 

As a basis of interpretation of the spectra ob- 
served with AKARI, we determine basic physi- 
cal parameters of six well observed objects. For 
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this purpose, we apply model photospheres as 
in the analyses of ordinary stars. However, the 
model photospheres of ultracool dwarfs including 
brown dwarfs are by no means well established 
yet. It has been discussed that dust grains in 
the atmosphere play an important role in the for- 
mation of spectra of ultracool dwarfs: especially 
spectra of L-type dwarfs (T off = 2200-1400 K) 
are strongly influenced by the continuous opac- 
ity of dust, while spectra of even cooler T-type 
dwarfs are dominated by deep molecular absorp- 
tion bands (|Tsuii et al. Ill996l ). In particular, we 



have no definite guideline for treating dust forma- 
tion in the photosphere. Some attempts to model 
the photospheres of brown dwarfs with dust clouds 
have been carried out by several groups and t hese 
are discussed in detail bv lHelling et al.l ( 2008f ). It 
was shown that all models agree on the global 
structure of the dusty photospheres, but predicted 
features such as the emergent spectra differ con- 
siderably, depending on the opacity data and re- 
flecting the extreme complexity of such objects. 

We apply the Unified Cloudy Model (UCM: 
Tsuiil 120021 l2005h as an example of such dusty 
model photospheres to the analysis of our objects. 
We assume Solar metallicity throughout, and es- 
timate T e ff, T cr (critical temperature that defines 
the upper boundary of the dust cloud; see Sec- 
tion 14.11 for detail), and log g from the relative 
SEDs (i.e. shapes of the spectra). We then fit 
the observed and predicted spectra on an absolute 
scale and estimate the radii of individual objects 
with the use of the known parallaxes. Based on 
the final fits of the model spectra to the observed 
ones, we try to interpret the new features in the 
spectral region between 2.5 and 5.0 /xm observed 
with AKARI for the first time. 

2. Observations and Data Reduction 

The AKARI /IRC covered its wavelength range 
with three independent cameras operated simulta- 
neously, namely the NIR (near-infrared), MIR-S 
(mid-infrared short), and MIR-L (mid- infrared 
long) channels. Our observations were all car- 
ried out in the AOT (Astronomical Observa- 
tion Template) IRC0 4 with parameters of "b;Np" 



R = A/AA = 120 at 3.6 /jm (jOhvama et al.ll2007l) . 

Targets were placed in the 1' x 1' point source 
aperture mask to minimize the contamination 
from nearby stars and background sky. The aper- 
ture mask was only available for the NIR channel, 
and the spectra obtained in this programme were 
limited to the near-infrared wavelength range of 
2.5-5.0 ^m. 

Our target list consisted of 30 brown dwarfs 
selected by their expected fluxes (to be bright 
enough for the AKARI /IRC instrument to pro- 
vide high-quality spectra within reasonable num- 
ber of pointings) and their spectral types (to sam- 
ple various types from L to T). A pointed observa- 
tion by AKARI allowed about 10 minutes of ex- 
posure. At least two observations were requested 
per target to obtain data redundancy. Due to se- 
vere visibility constraint of the satellite, not all 
requested observations were carried out. There 
were 18 pointing opportunities for 11 objects in 
the cryogenic mission phase. A summary of the 
observed targets is given in Table [IJ and the obser- 
vation record is presented in Table [2] The tables 
include observations which were not successful, for 
completeness. 

The standard software toolkit IRCSPEC.RED 



(Ohv ama et al.ll2007l) was used for the data reduc- 
tion^. The processing was performed in January 
2008 with the pre- released version of the toolkit, 
which should be equivalent with the public version 
"20080528" . We followed the standard reduction 
recipe. Fine tuning of the on-source / off-source 
mask on the spectral images was adopted to obtain 
the maximum signal and minimum contamination. 
Corrections of instrumental effects and wavelength 
/ flux calibrations were all done automatically in 
the toolkit, and we did not observe any calibration 
sources by ourselves. Since our observations were 
in principle slitless spectroscopy, the major source 
of the wavelength error was the determination of 
the reference point. The typical error was 0.5 
pixel of the detector or ~ 0.005 /jm (| Ohvama et al 



(jLorente et al.l 120081 ). In this mode a grism was 



20071 ). but could be larger on some occasions. We 
applied a small correction (—0.008 fim) to 2MASS 
J152322+3014 by comparing the position of CH4 
Q-branch feature with other objects and with the 
synthesized spectrum. The overall flux calibration 



used to achieve the highest available spectral res- 
olution; the dispersion was 0.0097 /im/pixel or 



lr The software is available from the AKARI observers' site; 
http: //www. ir . isas . jaxa . jp/AKARI/Observation/ 
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Table 1 

Summary of the targets in the current programme 



Object Name 


R.A.(J2000) 


Dec.(J2000) 


Sp.Type 


[J] 


[H] 


[K] 


W] 


ref. 


2MASS J04151954-0935066 


01 


15 


19 


54 


-09:35:06.6 


T8 


15 


32 


15 


70 


15 


83 


13.28 


1.3 


SDSS J053951. 99-005902.0 


05 


39 


52 


00 


-00:59:01.9 


L5 


13 


85 


13 


01 


12 


10 


11.32 


1.3 


2MASS J05591914- 1404488 


05 


59 


19 


14 


-14:04:48.8 


T4.5 


13 


57 


13 


64 


13 


73 


12.14 


1.3 


SDSS J083008. 12+482847.4 


08 


30 


08 


25 


+48:28:48.2 


L9 


15 


22 


14 


40 


13 


68 


11.98 


1.3 


2MASS J12171110-0311131 


12 


17 


11 


10 


-03:11:13.1 


T7.5 


15 


56 


15 


98 


15 


92 


13.96 


1.3 


GD 165B 


14 


24 


39 


09 


+09:17:10.4 


L3 


15 


64 


11 


75 


11 


09 


12.93 


1.3 


SDSS J144600.60+002452.0 


14 


46 


00 


61 


+00:24:51.9 


L5 


15 


56 


14 


59 


13 


80 


12.54 a 


1.3 


2MASS J15232263+3014562 


15 


23 


22 


63 


+30:14:56.2 


L8 


15 


95 


15 


05 


11 


35 


12.86 


1.5 


2MASS J17114573+2232044 


17 


11 


45 


73 


+22:32:04.4 


L6.5 


17 


09 


15 


80 


11 


73 


N/A 


2 


SDSS J175032.96+175903.9 


17 


50 


32 


93 


+17:59:04.2 


T3.5 


16 


14 


15 


94 


16 


02 


14.95 b 


1,4 


e Ind Ba+Bb 


22 


04 


10 


52 


-56:46:57.7 


T1+T6 


11 


91 


11 


31 


11 


21 


9.97 b 


2.1 



Refer ences. - [J, H K] (11 iKnapp et al.l |[2004D . 

[L'\ f4) lPatten et al.l ||2006D . (5) iLeggett et al.l [|2002h 

a Estimated magnitude. 

h Spitzer /IRAC 3.5 /im band magnitude. 



(2) 2MASS Point Source Catalog, (3) Golimowski et al ] l|2004h 



error is 10 per cent in the middle of the wavelength 
rang e and 20 per cent at the short / long edges. 
See lOhvama et al.l (|2007l ) for more details about 
data reduction and calibration. 

3. Near-Infrared Spectra of Brown Dwarfs 

Figure [T] shows the observed AKARI spectra 
of brown dwarfs in the sequence of their spectral 
types from L (bottom) to T (top). When an object 
was observed twice, the two data sets were pro- 
cessed independently and are plotted in different 
colors. In general, the two spectra are consistent 
with each other, demonstrating the high quality 
and reliability of the data. 

The overall SEDs of the observed spectra agree 
well with the predictions from the previous pho- 
tometry and theoretical works. The spectra of 
L-type dwarfs are rather smooth and featureless, 
mildly peaked at 3.8 /Ltm as a result of broad ab- 
sorption bands of H 2 at shorter wavelengths and 
CO at longer wavelengths. On the other hand, 
spectra of T-type dwarfs are dominated by deep 
molecular absorption features. As expected, we 
detect CH4 and H2O molecular bands. The bands 
deepen as the spectral type evolves. The band 
shapes are heavily distorted due to saturation in 
the latest object, 2MASS J041519-0935 (T8). 

The most remarkable finding in the current 
AKARI spectra is a detection of the CO2 fun- 
damental stretching-mode band at 4.2 [im in the 
late-L and T-type dwarfs. The identification of 
this band is demonstrated in Figure [2j in which 
the spectrum of 2MASS J055919-1404 is com- 
pared with synthesized spectra of C02- Spectra 
of CO and PH3 are also presented for reference. 
The shape of the band profile, especially a sharp 
edge around 4.2 confirms the identification of 
the feature as CO2. 

This is the first detection of CO2 in brown 
dwarfs, since the wavelength is inaccessible from 
the ground or even from airborne observations due 
to huge opacity of the molecule itself in the telluric 
atmosphere. The band is most apparent in the 
T4.5 dwarf 2MASS J055919-1404, and is clearly 
recognized in 2MASS J041519-0935 (T8) and e 
Ind Ba+Bb (T1+T6). It is seen in the L9 dwarf 
SDSS J083008+4828 and very marginally in the 
L8 dwarf 2MASS J152322+3014. We suspect that 
the formation of the molecule takes place in a rel- 



atively low temperature environment. 

15 




3.5 4.0 4.5 5.0 

Wavelength [/xm] 



Fig. 2.— The spectrum of 2MASS J055919-1404 
is compared with model spectra of CO2, PH3, and 
CO. The model spectra are calculated by a simple 
plane-parallel configuration with column densities 
and excitation temperatures of; 10 18 cm -2 and 700 
K for C0 2 , 10 19 cnT 2 and 700 K for PH 3 , and 10 18 
cm~ 2 and 1000 K for CO, respectively. The C0 2 
and CO data are taken from HITEMP database 



(jRothman et all 119971). and PHq data from HI- 
TRAN2008 (jRothman et al.l l2009h . The model 
spectra are normalized, and scaled and shifted for 
comparison. 

Another important result from our observations 
is that we see the CO fundamental band in the 
4.4-5.0 /Ltm region in all of our samples including 
the coolest T8 dwarf, 2MASS J041519-0935. The 
detection of the CO mo lecule in a late-T dwarf 
was fi rst re ported by Noll. Geballe. fc Marlev 
(|l997l ) and lOppenheimer et al.l (|l998h in the 
T6 dwarf Gl 229B. Tw o more objec t s hav e 



been added recently by Geballe et al. ( 20091 ) 



2MASS J09373487+2931409 (T6) and Gl 570D 
(T7.5). The presence of the CO band is also 
sug gested from the Spitzer 3.5-7.9 /xm pho tome- 
try (|Patten et al.ll2006HLeggett et al.ll2007t ). Our 
AKARI spectra provide direct and indubitable 
evidence of the molecule and strongly support the 
idea that the presence of the CO band is a com- 
mon feature in the latc-T dwarfs. 
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Table 2: Observation log 



Object Name 


Date 




OBSID 


Remarks 


2MASS 7041 519— 0935 


200 (■ 


(Jz- 


18 


1 toaaat 

1 ( ZUUUo- 


AA1 

UU1 


Crlicmtin v 

VJ 11 WO UlllfcL 


2MASS J041519-0935 


200 (■ 


(Jz- 


18 


1 ( ZUUUo- 


AAO 

UUz 


frhnstiTier 


2MASS J041519-0935 


2007- 


08- 


23 


CT 1 O C AO A 

5125080- 


AA1 

001 


Rc- observation 


2MASS 1041519—0935 


2007- 


08- 


24 


ri O r A O 1 

5125081- 


001 


R e-nbservn tion 


SDSS T053951-0059 


2006- 


A A 

09- 


17 


1 7AAAAA 

1720009- 


A A 1 

001 




2MASS T055Q1Q— 1404 


zuuu- 




22 


1 79000fi 






2MASS J055919-1404 


2006- 


09- 


22 


1720008- 


001 




SDSS J083008+4828 


2006- 


10- 


20 


1720007- 


001 




SDSS J083008+4828 


2006- 


10- 


21 


1720007- 


002 




2MASS J121711-0311 


2007- 


06- 


26 


1720068- 


001 


Too faint 


GD165B 


2007- 


07- 


24 


1720074- 


001 


Data lost 


SDSS J144600+0024 


2007- 


08- 


02 


1720072- 


001 




eps Ind Ba+Bb 


2006- 


11- 


02 


1720003- 


001 




eps Ind Ba+Bb 


2006- 


11- 


02 


1720004- 


001 




2MASS J152322+3014 


2007- 


01- 


26 


1720002- 


001 




2MASS J171145+2232 


2007- 


03- 


05 


1720001- 


001 


Too faint 


SDSS J175032+1759 


2007- 


03- 


17 


1720050- 


001 


Confusion 


SDSS J175032+1759 


2007- 


03- 


17 


1720050- 


002 


Confusion 



We also detect the CH4 molecule in SDSS 
J053951 — 0059 (L5 ) , co nfirming the previous re- 
port by iNoll et al. I (l2000h that the molecule is al- 
ready present in the atmosphere of L5 dwarfs. The 
feature is not clear in another L5 source SDSS 
J144600+0024, partly because of its relatively low 
S/N. The band develops in the L8 dwarf, 2MASS 
J152322+3014, and starts showing the P- and R- 
branches. 

Because of moderate spectral resolution it is dif- 
ficult to identify any more molecules in the spec- 
tra. 

3.1. Comments on individual objects 

3.1.1. SDSS J 05 3 951 - 0059 



Fan et al. (|2000f) selected this object as a high- 



z quasar candidate from the SDSS commissioning 
image and found it to be an L5 dwarf through 
follow-up spectroscopy. Due to its relatively high 
flux level (~ 8 mJy at 3.8 /im) the quality of the 
obtained spectrum is good. We can clearly recog- 
nize the presence of CH4 Q-branch absorption at 
3.3 fim on the otherwise rather smooth spectrum. 



3.1.2. SDSS J 144 600+ 0024 

This L5 object was no minated as a brown dwarf 
by iGeballe et al. ( 2002 ). Because of relatively 
low S/N compared to another L5 object SDSS 
J053951— 0059 above, we do not recognize the CH4 
Q-branch dip in the spectrum of this source, while 
the overall SEDs are consistent with each other. 

3.1.3. 2MASS 3152322+3014 



McLean et al. (l2000h found that this source is 
of the latest L-type an d assigned it as L8/ L9. It 



was confirmed as L8 bv lGeballe etaD (l2002h . This 
source, together with SDSS J144600+0024, is one 
of the faintest objects among our current dataset, 
and the quality of the spectrum is not excellent. 
However, we clearly see a more developed CH4 
band than in the L5 sources. The CO2 band is 
marginally detected. 

3.1.4. SDSS .1083008+4828 



Geballe et al 



(|2002l ) classified this object as 
L8-T0, or L9 as an average. It is on the bor- 
der of L/T spectral types and is a good example 
to investigate the transition from L to T. There 
are rather remarkable changes in the spectral fea- 
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C02 



"e Ind Ba+Bb (T1+T6) 

CH4 




2.5 3.0 3.5 4.0 4.5 5.0 
Wavelength [pim] 



33 
c 



33 

p 

x 



33 

x 



Fig. 1. — The NIR spectra of brown dwarfs obtained by the AKARI /IRC. The spectra are ordered in 
the sequence of their spectral types from bottom (L5) to top (T8). When two observations were made 
for an object, the data were processed independently and the second spectrum is indicated in blue. The 
difference between the two observations represents the practical errors. Positions of major molecular bands 
are indicated. 
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tures between this source from the previous L8 
source 2MASS J152322+3014. We see that H 2 
and CO absorption bands become more prominent 
than earlier L-type sources. The CH4 Q-branch is 
very clear, but P- and i?-branches are less obvious. 
The CO2 band is clearly seen. 

3.1.5. elndBa+Bb 



This object was identified by IScholz etal 



(|2003l ) as a companion to the K5V star, e Ind A. 
Because of its proximity to the Sun (3.6 pc), 
it is the brightest among our observed targets. 
We observed the object twice within a day but 
found that one observation (OBSID=1720003- 
1) was significantly degraded by charged parti- 
cle hits before and during the exposure. The 
other spectrum (OBSID=1720004-1) shown in 
Figure [TJ clearly exhibits the four major molecular 
bands; CH4, H 2 0, CO and CO2. This object is in 
fact a binary system of Tl and T6 dwarfs, with 
if-band magnitudes of 14 and 16, respectively 
(|McCaughrean et all \2004) . AKARPs spectrum 
is a mixture of the two spectral types and we will 
not use it for quantitative analysis. 

3.1.6. 2MASS J055919-1404 

The object was identifi ed as a 'warm' T-dwarf 
by iBurgasser et al. ( 2000f) . From a near-infrared 
spectrum in 0.9-2.3 /mi they concluded that 
the source was near th e L/T transition border. 



Burgasser et al 



(|2006bh classified the object as 
T4.5. The AKARI spectra of this dwarf show very 
deep CO and CO2 absorptions beyond 4.0 /mi, to- 
gether with CH4 and H2O bands below 3.8 /jm. 
As a result the spectra show a peak at 4 /mi. A 
sharp peak, probably the CO band center appears 
very clearly between the R- and P-branches. Our 
ob servations are consis tent with the photometry 
bv lLeggett et al. I (|2002l) that the M' band flux of 
this source is fainter by a factor of three than that 
expected if CO/CH4 abundance follows thermal 
equilibrium. Our spectra confirm that CO is the 
major source of opacity in the wavelength region. 

3.1.7. 2MA SS J 041 51 9- 0935 

This object w as first identified as a T-type 



Burgasser et al. (|2006bh . It is one of the lat- 
est spectral-type brown dwarf, and is the coolest 
object in our sample. The first observations of 
this target with AKARI were carried out in 2007 
February. However, the data were considerably 
contaminated by ghosting from a nearby bright 
star. We reobserved the object a half year later in 
2007 August and obtained a clean spectrum. 

The spectra exhibit deep H2O and CH4 absorp- 
tion bands, which are almost completely saturated 
at the wavelengths below 3.5 /mi. On the other 
hand, the object has sufficient flux in the longer 
wavelength range, where we clearly identify the 
CO and CO2 bands. Two independent observa- 
tions show almost identical shape in this wave- 
length range, strengthening the reality of the de- 
tections. 

4. Comparisons of the Observed and Model 
Spectra 

In this section, we examine to what extent 
the newly observed spectra of the 2.5-5.0 /mi re- 
gion can be explained with the model spectra 
and what remains unexplained with the present 
models of the photospheres of cool dwarfs. For 
this purpose, we apply our Unified Cloudy Model 
(UCM), a brief description of which is given in 
Section 14.11 We assume local thermodynamical 
equilibrium (LTE) throughout this section. Then 
we outline our method of analysis with some exam- 
ples in Section 14.21 and the results of our analysis 
in Section [4731 The physical parameters estimated 
by our analysis are summarized in Section f4.4l Fi- 
nally, we examine the errors and limitation of our 
analysis based on the UCM in Section 14751 



brown dwarf by IBurgasser et al 



(|2002h based 
on near infrared spectroscopy. The object has 
been designated a "standard" T8 source by 



4.1. Predicted spectra of dusty dwarfs 

Generally, stellar spectra can be interpreted in 
terms of T e g, log g, chemical composition, and 
micro-turbulent velocity. However, a new feature 
in the spectra of brown dwarfs, compared with 
usual stellar spectra, is the effect of dust clouds 
formed in the photosphere. If the properties of 
dust clouds can be defined uniquely from the four 
basic parameters considered above, the spectra 
of brown dwarfs could eventually be interpreted 
in terms of these four parameters, namely T e ff, 
log 5, chemical composition, and micro-turbulent 
velocity. At present, the dependence of the dust 
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cloud properties on the four basic parameters is 
unknown, because the formation and disappear- 
ance of dust clouds in the photospheres of brown 
dwarfs is not yet well understood. 




Fig. 3. — The predicted spectrum of a UCM with 
T cff = 1300 K, T cr =1800 K, logg = 4.5, £ micro =1 
kms -1 , and Solar metallicity, showing the effect of 
the CO2 4.2 /Ltm band based on the line list (solid 
line) and on the band model opacity (dotted line) 
compared with the case without CO2 lines (dashed 
line) . 

On the other hand, some observations indicate 
that the properties of the dust clouds do not nec- 
essarily depend uniquely on such basic parame- 
ter as T e g. For example, infrared colors such as 
J — K plotted ag ainst T B ff show a large variatio n 
at a fixed T off (e.g. iMarlev et al.ll2005t rrsulill2005l) . 
Since these infrared colors depend directly on the 
properties of the dust clouds, this fact implies that 
the properties of the dust clouds should be differ- 
ent even for the same T e g. Also, one important 
observational result is that T e g shows little change 
over the spectral types between about L5 and 



T5 ( e.g. iGolimowski et al.l l2004t iNakaiima et al 



2004). In other words, the characteristics of the 
spectra used as signatures of the spectral types 
arc quite different even for the same T c g and 
hence they are not necessarily be determined by 
T c ff. This result implies that the basic features of 
the spectra of brown dwarfs should be determined 
by some additional parameter(s) other than T e ff. 
Again dust should play a major role in this respect 



as in the case of the infrared colors noted above. 

Thus, a major problem is how to consider the 
effect of dust in the interpretation and analysis of 
the spectra of brown dwarfs. For this purpose, we 
apply our Unified Cloudy Model (UCM) in which 
the dust forms at a layer where temperature is 
equal to the condensation temperature, T con d, but 
disappears at somewhat lower temperature which 
we refer to as the critical temperature, T cr . At 
this temperature, the radius of the dust grain r gr 
reaches the critical radius r cr , where the Gibbs en- 
ergy of condensation is maximum. Then, at this 
point, dust starts to grow larger and larger, and 
will precipitate from the gaseous mixture. For this 
reason, dust disappears in the layers with T < T cr , 
and will exist only in the layers of T cr < T < T cont x- 
Thus dust is present in the layers whose photo- 
spheric temperatures are within the above range, 
and this means that a thin dust cloud forms in 
the photosphere of bro wn dwarfs, as dis cussed in 
more detail elsewhere (|Tsuii|[2001 . 2002). It is to 
be noted that T con d is defined by the thermody- 
namical data and is hence full y controlled b y T e s 
and log g (see e.g. Figure 3 in iTsuiil (|2002t )'l. On 
the other hand, T cr is not predictable by any phys- 
ical theory at present. Instead, T cr is introduced 
as a free parameter in our UCM, and it serves as 
a measure of the thickness of the dust cloud (the 
dust cloud is thicker if the deviation of T cr from 
T cond is larger). 

With the introduction of T cr , the UCM is char- 
acterized by the five parameters including T cr in 
addition to the usual four parameters, T c g, log g, 
chemical composition, and micro-turbulent veloc- 
ity. From observations, we showed that at least 
one additional parameter is needed to describe the 
observed characteristics of brown dwarfs, and we 
now propose to identify T cr as the additional pa- 
rameter required from observations. 

A new feature in the AKARI observations dis- 
cussed in Section [3] is the detection of the CO2 
band at 4.2 ^m. Although CO2 is included in the 
construction of the UCM as well as in our eval- 
uation of the spectra with the use of the b and 
model opacity (e.g. Appendix of ITsuiil 120021 ). we 
recompute the spectra w ith the CO2 line list fr om 



the HITEMP database ( Rothman et al.l 119971 ) in 



view of the detection of the CO2 band in several 
objects. As an example, a predicted spectrum in- 
cluding CO2 based on the line list (solid line) is 
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compared with that without CO2 (clashed line) in 
FigureElfor a UCM with T cff = 1300 K, T cr = 1800 
K, logg = 4.5, £ m icro = 1 kms -1 , and Solar abun- 
dance. For comparison, the result based on the 
band model opacity is shown by the dotted line. 
We use the band model opacity during the itera- 
tions of model construction but apply the line list 
in all the computation of the spectra. So far, CO2 
was largely neglected in spectroscopy and model 
photospheres of cool dwarfs, but the result shown 
in Figure [3] reveals that the effect of CO2 absorp- 
tion is significant even within the framework of the 
LTE analysis. Since carbon atoms are mostly con- 
sumed in CO and/or CH4 in cool dwarfs, CO2 can- 
not be very abundant, but rather large /-values of 
CO2 transitions make the molecule an important 
absorber in the spectra of cool dwarfs. 

The major spectral feature in the 2.5-5.0 /zm 
region is the CH4 fundamental band at 3.3 fim, 
for which we apply the line list by Freedman 
(2005, private communication) discussed in de- 
tail by [fteedman^arieylALodderi |200i). Cer- 
tainly the CH4 line lists including the high exci- 
tation l ines are not yet satisfactory as discussed 
before (jTsuiil 120051 ) , but we believe that Freed- 
man's list is the best one currently available for 
the fundame ntal band. Other line l i sts in cluded 
are: H-,0 (IPartridge fc Schwenkel Il997h. CO 



( Guelachivili et all 19831: IChackerian fc Tipping! 
19831). OH (Ijacauinet-Husson et all Il999l). CN 



(ICernv et allll978UBauschlicher et alJll988l). and 



SiO flLavas. Maki fc 01sonll98ll:lTipping fc Chackerian 
198ll ). Also. NH^. PFFi. and H9S are considered on 
the basis of the band models. With these spectro- 
scopic data, the spectra between 2.5 and 5.0 /xm 
for our UCMs are evaluated at a resolution of 
0.01 cm -1 and convolved with the slit function of 
FWHM = 3000kms _1 , which is about the resolv- 
ing power of the AKARI spectrometer. Note that 
the observed spectra shown in Figures IH [5j & |6] 
are smoothed with R = 100, and that the error 
weighted average is taken if two valid spectra are 
available. 

4.2. Method of analysis and some exam- 
ples 

We restrict our analysis to the spectral region 
between 2.5 and 5.0 fim, which we have observed 
completely for the first time. Although analyses 
of different spectral regions do not necessarily pro- 



vide the same answer as to what are the best pa- 
ra meters of the photosp here, as shown in detail 
by ICushing et al. ( 20081 ). we restrict our analysis 
to the 2.5-5.0 /im region for this very reason. In 
fact, for our purpose stated at the beginning of 
this Section, it is meaningless to find a solution 
that provides a good fit to other regions but not 
to the 2.5-5.0 /itm region. We return to this prob- 
lem in Section 1431 

We assume that the object s we have observed 
are a l l of Solar metallicity ( Anders &: Grevesse 
19891 : lAllende Prieto et alJ 120021 ) and that the 
micro-turbulent velocity is near solar (1 kms -1 ) 
throughout. Then, major parameters that define 
the characteristics of a spectrum are T e ff, T cr , and 
logg. In fitting the observed spectrum with the 
predicted one, we refer to the empirical data sum- 
m arized in Ta b le El We start with the parallaxes 
bv lVrba et all (I200I which covered our full sam- 
ple. We notice that their values agree quite well 
with other results in general. We choose T e ff to be 
Ti g close to the so-ca lled empirical T e g suggested 
bv IVrba et all (|2004l ) (Tabled and examine the 
cases of T eS = T e ° ff and T e ° ff ± 100 K. As a result, 
we examine three values of T e g near the empirical 
T e g. But this choice of trial T e ff's is not neces- 
sarily successful as described below (e.g. a case 
of 2MASS J055919-1404 to be discussed in Sec- 
tion 14. 2p . This fact suggests that the so-called 
empirical T G ff cannot be fully realistic as will be 
discussed in Section 14.4.11 We also consider three 
cases of T cr = 1700, 1800, & 1900 K and three 
values of logg = 4.5, 5.0, 5.5. Thus, we consider 
3 x 3 x 3 = 27 combinations of these parameters 
(T e ff , T cr , and logg), and select the case that shows 
a best fit to the overall shape of the 2.5-5.0 /Ltm 
SED. 

We discuss the procedure outlined above for the 
case of SDSS J083008+4828 as an example. The 
empirical T cS of this L9 dwarf is 1327 K (Tabled, 
but the CH4 3.3 /jm band is too weak for such 
a low value of T or f. Therefore we start our trials 
with T c fr = 1300 K and extend then to higher T e g, 
rather than to examine the cases of T e g = 1300 ± 
100 K. It is rather difficult to find a good fit in 
this object, but one possible solution is obtained 
for (T off , T cr , logg) = (1500, 1700, 4.5) from our 
27 trial models. The effect of T c ff along with the 
possible best case is shown in Figure 2^,: only the 
case of T c ff =1500K provides a reasonable fit to the 
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P- and i?-branches of CH4 fundamentals, although 
the predicted Q-branch is a bit stronger than that 
in the observed spectrum. 

The effect of T cr at fixed values of T cS = 1500 
K and \ogg = 4.5 is illustrated in Figure [4}d. Al- 
though the fit to the observed CO2 4.2 jxm feature 
can be improved by higher values of T cr (i.e. with 
thinner dust cloud), the fit to the CH4 3.3 fim fea- 
ture appears to be worse. The effect of log g is 
rather modest as shown in Figure Ht, but the low- 
est gravity of log<? = 4.5 provides better fits both 
to the CH4 3.3 /im and CO2 4.2 /im bands than the 
higher values. The final fit we find, however, is not 
so good around 2.7, 4.2, and 4.5 /im and we will 
discuss these wavelength regions in Section |4~3"1 

Another example is the case of 2MASS J055919-1404. 
The empirical T cS of this T4.5 dwarf is 1469 K 
(Table[3]), and we, at first, examined the cases of 
T cff = 1400 ± 100 K. It was difficult to find a 
good fit, but we found that the case of T e g = 1300 
provides a relatively close fit to some features. 
However, an anonymous re feree called our atten- 
tion to a recent result by ICushing et al. I (l2008h 
who obtained T c s = 1200K for this object, and 
suggested to examine this case. Although this 
case differs by a large value (« 300 K!) from the 
empirical value, we now include T e g = 1200 K in 
our trial and find that this case actually provides 
a better fit as shown in Figure [S£i. For example, 
the relatively strong CH4 fundamental band can 
be better explained by T e g = 1200 K than by T e g 
= 1300 K. However, the CO2 4.2 /zm feature as 
well as the strong absorption at the positions of 
CO fundamentals cannot be fitted at all. We will 
return to this problem in Section 14.31 

If we assume values of T cr lower than 1900 K, 
the fits to the CH4 3.3 /jm band tend to become 
worse while those to the CO2 and possible CO 
bands show little change as seen in Figure [SJ). The 
J — K color near zero (Table [3]) also suggests that 
the dust cloud is rather thin, consistent with the 
high value of T cr . The effect of log g is examined 
for the fixed values of T eS = 1200 K and T cr = 
1900 K. It is found that the spectrum shows little 
change for the change of log g, but the 2.8-3.0 fj,m 
feature tends to fit better for the lower values of 
log g (Figure [5b) . We conclude that the case of 
(T cf f, T cr , \ogg) = (1200, 1900, 4.5) is the best 
compromise for this object. 



4.3. Results 

We now extend the same procedure outlined 
in Section 14.21 to the other objects except for e 
Ind Ba+Bb because of the composite nature of its 
spectrum due to its binarity. The observed and 
predicted spectra can be fitted rather well in the 
remaining four objects (see Figure H]), and thus 
they are more or less easier to analyze compared to 
the cases discussed in Section |4"T2"1 For this reason, 
we only give the resulting best fits for their^. The 
parameters that give the best fits for all six objects 
are summarized in Table |4] After we find the pos- 
sible best model for each object, we measure the 
vertical shift in logarithmic scale between the ob- 
served spectrum and the predicted one based on 
the best model we found. The emergent flux from 
the unit surface area of the object, F v (in units of 
ergcm~ 2 s _1 Hz -1 ), is written as 

\ogF v = log /„ - 2 \og{R/d) - 23.497, 

where /„ (in unit of Jy) is the observed flux, d 
is the distance to the object based on the mea- 
sured parallax, and R is the radius of the object. 
The vertical shift provides R/d. With the known 
distance, we can estimate the radius R and the 
result is also given in Table [4] The resulting fits 
in absolute scale are given in Figure |6] 

4.3.1. SDSS J 05 3 951 - 0059 

The empirical T e g of this L5 dwarf is 1690 K 
(Table [3]), and we find that the case of (T e g, T cr , 
log g) = (1800, 1800, 5.5) shows the best fit for 
this object. In fact, the overall shape of the SED 
as well as the strengths of the CH4 Q-branch at 
3.3 fim and the CO fundamental at 4.6 /im is rea- 
sonably fit for this case as shown in Figure [6^,. It 
is to be noted that the effect of T cr is of compara- 
ble importance to that of T e g. In fact, the column 
densities of both dust and molecules generally in- 
crease at the lower T e g (at fixed T cr ) and also at 
the lower T cr (at fixed T c g). For this reason, the 
effects of T e g and T cr on the spectra are often dif- 
ficult to discriminate. We have considered the ef- 
fects of both T e ff and T cr from the beginning, and 
looked for the best solution. The effect of log g on 



2 The 
tra 



numerical details of the predicted spec- 
and the UCMs are available from 
http: //www.mtk. ioa. s .u-tokyo . ac . jp/~ ttsuj i/export/ucmLM, ucm 
respectively. 
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■ SDSS J083008+4828 (L9) 

■ a) Tcr=1700 




I i i I i i i i I i i i i I i i I 

3 4 5 

A(/um) 

Fig. 4. — The observed spectrum of SDSS J083008+4828 (L9) is compared with the predicted spectra based 
on the UCM. The best fit obtained for (T off , T cr , logg) = (1500, 1700, 4.5) is shown in each plot. Flux is 
adjusted by shifting the model spectra vertically. The same amount of shift is applied to the other predicted 
spectra. Hence no attempt is made to fit them to the observed spectrum, but are shown here for comparison 
purposes, a) Effect of T c g under the fixed values of T cr =1700K and \ogg — 4.5. b) Effect of T cr under the 
fixed values of T e g=1500K and logg = 4.5. c) Effect of logg under the fixed value of T cr =1700K and T cff = 
1500K. 
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1 1 1 1 1 1 1 — 

- 2MASS J055919-1404 (T4.5) 




j i i i i i 

3 4 

A(/im) 



Fig. 5. — The observed spectrum of 2MASS J055919— 1404 (T4.5) is compared with predicted spectra based 
on the UCM. The best fit is obtained for (T cff , T cr , logg) = (1200, 1900, 4.5). See the caption of FigureH 
as for other details, a) Effect of T c g under the fixed values of T cr =1900K and log g = 4.5. b) Effect of T cr 
under the fixed values of T e ff=1200 K and log g = 4.5. c) Effect of \ogg under the fixed value of T cr =1900 
K and T cS = 1200 K. 
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Table 3: The known observed data 



no. 


object 


sp. type a 


J - K b 


parallax c 












Opt./IR 




(marcsec) 




(K) 




1 


SDSS J053951-0059 


L5/L5 


1.45 


76.12 ±2.17 


1690 (- 


hl73, 


-137) 


2 


SDSS J144600+0024 


L6/L5 


1.76 


45.46 ±3.25 


1592 (- 


1-175, 


-140) 


3 


2MASS J152322+3014 


L8/L8 


1.60 


57.30 ±3.27 


1330 (- 


1-142, 


-112) 


4 


SDSS J083008+4828 


L8/L9 


1.54 


76.42 ±3.43 


1327 (- 


1-140, 


-111) 


5 


2MASS J055919-1404 


T4.5 


-0.16 


95.53 ± 1.44 


1469 (- 


1-153, 


-122) 


G 


2MASS J041519-0935 


T8 


-0.51 


174.34 ±2.76 


764 ( 


+88, 


-71) 



Ref 
a) 
b) 



■L. 

Kirkpatrickl ()2005l) for L and Burgassc r et alj (I2006M for T 



Knapp et all (|2004h 



c) IVrba "et al.1 (l2004h 



Table 4: Basic parameters estimated from the model fittings using UCMs 



no. 


object 


T off (K) 


T CI (K) 


logs' 


R/Rf 


AT off (K) b 


1 


SDSS J053952-0059 


1800 


1800 


5.5 


0.804 


-110 


2 


SDSS J144600+0024 


1700 


1700 


4.5 


0.716 


-108 


3 


2MASS J152322+3014 


1500 


1700 


4.5 


0.684 


-170 


4 


SDSS J083008+4828 


1500 


1700 


4.5 


0.700 


-173 


5 


2MASS J055919-1404 


1200 


1900 


4.5 


1.175 


-269 


6 


2MASS J041519-0935 


800 




4.5 


0.767 


-36 



a) radius R relative to the Jupiter's radius Rj (see Section [ 

b) AT cff = T cff (Table [3] by empirical method) - T eS (Table H by our model fitting). 
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the spectrum is rather minor compared with the 
other parameters, but the lower gravities may be 
reasonably excluded from the overall poor fits. 

The resulting fit in absolute scale (Figure (S^i) 
shows that the spectrum observed by AKARI can 
reasonably be accounted for by our model. 

4.3.2. SDSS J144600+ 0024 

Although the empirical T e s of this L5 dwarf 
is near 1600 K (Table [3]), we find it very difficult 
to assume T e s < 1600 K since the observed CH 4 
feature at 3.3 /im is very weak if present. Again by 
the procedure outlined in Section |4~21 we find that 
a case of (T cff , T cr , log g) = (1700, 1700, 4.5) shows 
the best fit for this object. For example, we see 
that not only the strength of the CH4 3.3 /im band 
but also the overall shape of SED are best fitted 
with T cff =1700K. The effect of T cr is rather large 
and we can easily exclude such high values of T cr 
as 1800 and 1900K. Also, the very red J — K color 
(Table [3]) suggests that the dust cloud is pretty 
thick and T cr cannot be as high as 1900K. It is 
found that the overall SED is again not sensitive 
to log g, but the strengths of both the CH4 Q- 
branch at 3.3 /im and CO fundamental band could 
discriminate the best log g of 4.5 from the higher 
values of log g = 5.0 and 5.5. 

The resulting fit in absolute scale in Figure |6Jd 
shows that the observed spectrum of the 2.5- 
5.0 /xm region, which possesses almost no distinct 
features, can reasonably be accounted for by our 
model. 

4.3.3. 2MA SS Jl 52322+ 3014 

The empirical T c g of this L8 dwarf is 1330K 
(Tabled]), but we found that the best fit is found 
for (T off , T cr , logs) = (1500, 1700, 4.5). In fact, 
the modest strength of the CH4 3.3 /mi band sug- 
gests that T c ff cannot be below 1300 K, and it is 
only the case of T c s = 1500 K that provides a rea- 
sonable fit not only to the Q- but also to the P- 
and i?-branches of the CH4 fundamental band. We 
confirm that such a fit as shown for T cr = 1700 K 
can no longer be found for T cr as high as 1800-1900 
K. This result is consistent with the red J—K color 
of 1.60 (Table[3]), which implies that the dust cloud 
of this L8 dwarf is rather thick. The strength of 
the CO fundamental band is almost independent 
of logs, but the P- and i?-branches of the CH4 



fundamentals depend on log g and can only be ex- 
plained with the low gravity of logs = 4.5. The 
effect of logs 011 the infrared spectrum is again 
rather modest and accurate determination of log g 
is difficult. 

A detection of CO2 at 4.2 /im suggested in Sec- 
tion [3] can be well supported by our model predic- 
tion that shows a depression at the expected posi- 
tion of the CO2 4.2 /im band and even shows a rea- 
sonable fit quantitatively (see Figure [5b). The fit 
in the region of the CO fundamental band appears 
to be poor and the deep absorption features near 
4.5 and 4.9 /xm remain unexplained. We suspect, 
however, that the CO fundamental band provides 
the major contribution as in later type dwarfs to 
be discussed in Section [A. 3. 41 

Thus the spectrum observed by AKARI now 
shows some dissonances with the predicted spec- 
trum in this late-L dwarf as shown in Figure [6b. 

4.3.4. SDSS J083008+4828 

As already discussed in some detail in Sec- 
tion H2J we find {Tea, T cr , logs) = (1500, 1700, 
4.5) for this object. The final fit is shown in Fig- 
ure [6b! . We notice that the resulting parameters 
are the same as that for 2MASS J152322+3014 
discussed in Section 14.3.31 This result is also con- 
sistent with the similarity of the empirical T e s and 
of the J — K color between these two late-L dwarfs 
(see Table [3]). Nevertheless, the observed spectra 
of the two objects differ significantly. For exam- 
ple, the observed feature near the 2.7-2.9 /nm re- 
gion (largely due to H 2 ^1,^3 and 2v 2 bands) 
shows a large difference between these two objects 
of the same physical parameters. Since this feature 
is well reproduced by our models in other objects 
except for SDSS J083008+4828 (see Figure!]), this 
should not be due to any systematic effects in our 
models. 

A more interesting result is that the observed 
CO2 4.2 /im feature in this object is very strong as 
noted already (Section[3]), and that it is too strong 
compared with the prediction in marked contrast 
to the case of 2MASS J152322+3014. It appears 
that the CO2 4.2 /im feature cannot be explained 
by any combination of T e g, T cr , and logs with our 
UCM essentially based on the LTE assumption, 
while those of 2MASS J152322+3014 may be more 
consistent with the LTE prediction. Such a large 
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Fig. 6. — The observed and predicted spectra based on the UCMs are compared on an absolute scale: 
a) SDSS J053951-0059 (L5); b) SDSS J144600+0024 (L5); c) 2MASS J152322+3014 (L8); d) SDSS 
J083008+4828 (L9); e) 2MASS J055919-1404 (T4.5); and f) 2MASS J041519-0935 (T8). 
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difference in CO2 band in the two objects of sim- 
ilar physical parameters cannot be explained by 
differences in dust cloud properties if, as we find, 
T cr is nearly the same around 1700K in both ob- 
jects. The deep absorption at 4.5 fxm can be iden- 
tified with the i?-branch of the CO fundamental, 
but cannot be explained within the framework of 
the LTE assumption as well. 

Thus, there are serious difficulties in under- 
standing the spectrum observed by AKARI in 
that the region between 4 and 5 /jm cannot be 
fitted at all by the LTE prediction based on our 
UCM. 

4.3.5. 2MA SS J 05591 9-1404 

Details of the model fitting of this object were 
described in Section l4~2| and our final fit on an ab- 
solute scale is shown in Figure We are shocked 
to see that the discrepancy between the observed 
and predicted spectra is so large in the region be- 
yond 4.1 /j,m. The region of CO i?-branch is too 
deep to be explained by our LTE model and CO 
P-branch also appears to be appreciable. Thus we 
confirm that the CO bands are actually enhanced 
not only in late-T dwarfs but also in mid-T dwarfs, 
as anticipated by p hotometric observations (e.g. 
Leggett et al1l2007l i. 

In this object, a more important result is that 
the new spectrum observed by AKARI provides 
definite evidence for the presence of the strong 
CO2 band in brown dwarfs for the first time. The 
depression at 4.2 /im is too large to be accounted 
for by the predicted CO2 band and the discrep- 
ancy between the observed and predicted spectra 
is quite serious. All these new features due to CO 
and CO2 cannot be accounted for by the present 
LTE model such as our UCM. In contract to these 
serious disagreements in the region beyond 4 /im, 
the observed and predicted spectra agree quite 
well in the region between 2.5 and 4.1 jum. 

Recent analysi s of t he SED of this object 
by ICushing et all (|2008l ) concluded that T cS = 
1200 K, log g = 5.5 and / scd = 4, where / sod 
is the cloud sedimentation efficiency parame- 
ter and a larger value implies greater sedimen- 
tation efficiency resulting in a thinner cloud 
(jAckerman fe Marlevl 120011 ) . Thus / se d is a pa- 
rameter that may have a similar role as our T cr , 
with larger / se d corresponding to higher T cr . Thus 



their result showing relatively large value of / sc d 
■1 is consistent with our relatively high value of 
T cr = 1900 K, and we agree with them on the 
properties of the dust cloud in this T4.5 dwarf. 
We also agree that T e g = 1200 K for this object. 
Thus our results for the cloud properties as well 
as for r s ff show goo d agreements with those by 
Cushing et al.l (|2008l ). despite the fact that dif- 
ferent spectral regions are analyzed with different 
models. Unfortunately, the results for gravity dif- 
fer considerably. However, SEDs are not very 
sensitive to the gravity and this is the parameter 
most difficult to determine from the model fittings 
(as will further be discussed in Section 14. 4. 2|l . 

In summary, some dissonances between the ob- 
served and predicted spectra, seen in the late-L 
dwarfs, appear to be quite distinct in this mid-T 
dwarf. 

4.3.6. 2MASS J041519-0935 

The empirical T cS of this T8 dwarf is 764 K (Ta- 
ble [3} and, if a dust cloud forms, it is situated be- 
low the layer of the optical depth unity in such a 
very cool dwarf. For this reason, the spectrum no 
longer depends on T cr , as we have actually con- 
firmed. Accordingly, we apply the dust-free mod- 
els in which dust grains have all precipitated just 
after they are formed (this means T cr = T cond ), 
and we find the best possible fit for (T e ff, T CT , log g) 
= (800, T cond , 4.5). 

Determination of the physical properties of the 
coolest brown dwarfs including 2MASS J041519-0935 
has bee n discussed bv lBurgasser. Burrows fe Kirkpatrick 
(l2006aT) . who calibrated the strengths of H2O 
bands and relative fluxes in terms of T e g and log g. 
As a result, they suggested T e g = 740-760 K and 



=4.9-5.0 for this object. Also. lSaumon et al 



(|2007l ) suggested a similar result of T e g =725- 
775 K and \ogg =5.00-5.37 based on the analysis 
of a spectrum including new data obtained with 
Spitzer Space Telescope. Our result is reasonably 
consistent with these results, although the val- 
ues of log g differ from each other. Thus analyses 
of the different spectral regions by the different 
models show more or less consistent results in this 
object at least for T e g. 

The overall fit of the predicted spectrum of 
the possible best model is not so bad (see Fig- 
ure [6f), but we confirm that the region of the 
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CO bands cannot be fitted at all. Such a pos- 
sible large discrepancy in the predicted and ob- 
served CO bands in late-T dwarfs was previously 
suggested based on the ground based observa- 
tion of Gl 229B (iNoll. Geballe. fc Marlevl Il997t 
Qppenheimer et all 19981 ), and more recently fo r 
two other late-T dwarfs bv iGeballe et all (|2009l ). 
We confirm these previous results and add a new 
case of largely enhanced CO fundamentals in late- 
T dwarfs. We also note that the i?-branch is much 
stronger compared with the P-branch. 

A more interesting result is that the CO2 band 
at 4.2 /im appears quite distinctly in this late-T 
dwarf too. It is impossible to explain both the CO 
and also the CO2 bands by the LTE models, and 
we will discuss such cases in the Section [5] 

In conclusion, the new spectral region explored 
by AKARI for the first time reveals: first, the en- 
hancement of the CO fundamental bands is a gen- 
eral feature and, second, the CO2 band at 4.2 ^m 
appear to be an important absorber, in late-L and 
T dwarfs. 

4.4. Basic physical parameters of six brown 
dwarfs 

We now discuss briefly the resulting T e g, T cr , 
log 5, and R/Rj values summarized in Table 21 
Also, we discuss briefly the effects of other basic 
parameters, namely chemical composition (in Sec- 
tion I4.4.4|) and micro-turbulent velocity (in Sec- 
tion 033]), which we assumed to be the Solar. If 
we try to determine these two parameters by the 
same way as we have done in Section 14.21 and 14.31 
for the three basic parameters, we must consider 
at least 3x3x3x3x3 = 243 models instead 
of 27 models for each object, but this is clearly 
impractical. Moreover, we do not think that the 
chemical abundances and micro-turbulent velocity 
can be determined from such low resolution spec- 
tra as those we have at hand. 

4.. 4-1- Effective temperature and radius 

We notice first that our T e g's do not agree with 
the so-called empirical effective temperatures ob- 
tained from the observed lumin osities and a con - 
stant radius such as R = 0.9-RT ljVrba et al.ll2004l ) 
(see 7-th column of Table 0J. In fact, we al- 
ready notice through Section 14.3.11 and 14.3.61 that 
the observed SEDs cannot be fitted well with 



the predict e d one s assuming the empirical T e g by 



Vrba et al.l (J2004J). The reliability of the empirical 



T e ff largely depends on the accuracy of the values 
of the radius. The constant radius of R = 0.9 Rj 
was b ased on a statistical analysis by iBurgasser 



(200 If) of the evolutiona ry models of cool dwarfs 
by Burrows et al.l ( 19971 ). but it is unknown if it 
is appropriate to apply such a value to all the ob- 
jects. 

It is known that the radii of cool substellar ob- 
jects are independent of mass to within 30 per cent 
for a broad range o f masses from 0.3 to 70 Mj 
( Burrows et al-llioOll ). The resulting radii of all 
our six objects (Table [4]) appear to be consis- 
tent with the prediction of the evolutionary mod- 
els and confirm that the radii of ultracool dwarfs 
are within approximately 30 per cent of Jupiter's 
radius. The mean radius for our six objects is 
0.81 R j, and thus we agree with IBurgasser I (|200lj ) 
in that the radii of ultracool dwarfs are statisti- 
cally smaller than Jupiter's radius. However, the 
variation of radii among brown dwarfs appears to 
be appreciable (Table [4j and may not be repre- 
sented by a single mean value. For this reason, 
the so-called empirical effective temperature can- 
not be completely reliable. Although we use it as 
a guideline in our analysis (Sections 14.21 &: |4T3|) . it 
may even be misleading as we showed in the case 
of 2MASS J055919-1404 (Section [OJ . 

Also, because of such variations of the radii in 
real ultracool dwarfs, an attempt to analyze the 
observed spectra on an abso lute scale ba sed on an 
assumed radius of R = Rj ( Tsujil 20051 ) does not 
necessarily work well. For example, we showed 
before that the near infrared spectrum of 2MASS 
J152322+3014 observed with Subaru, reduced to 
an absolute scale with the known parallax and 
with the assumption of R = Rj, agreed with the 
predicted spectrum of T c g = 1300 K. For this rea- 
son, we sugg ested T e g of this object to be ~ 1300 K 
(|Tsuiill2005f ). However, we find that the shape of 
the AKARI spectrum can be best explained by 
T e g = 1500 K (Section I4.3.3|) . and that a consis- 
tency with the SED on an absolute scale can be 
achieved if we abandon the ad hoc assumption of 
R = Rj (see Figure Et). 

In conclusion, we believe that the observed 
spectrum (or relative SED) can provide a more 
direct and realistic estimation of the effective tem- 
perature of an individual object. Further, an em- 
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pirical estimation of the radius can be possible 
if the SED can be reduced to an absolute scale 
with the known parallax. With this conclusion, 
we now think it inappropriate to analyze the ob- 
served spectrum in absolute scale based on an ad 
hoc assump tion of R = Rj, as we proposed before 
(|Tsuiill2005l) . 



4-4- Gravity 

The determination of log g is more difficult since 
the low-resolution infrared spectra depend little 
on log g, as can be seen in our Figures 0] & El 
Also, this difficulty is clearly shown by a goodness- 
of-fit (statistic Gk defined by Equation (1) in 
Cushing et al. 1 (l2008h and shown in their Figure 5, 
in which Gk only shows minor changes for different 
logo's (possibly except for no cloud case). 

Inspection of Table|4]reveals that our result ap- 
pears to be biased towards low gravity, namely 
\ogg — 4.5 for 5 out of 6 objects. This may be 
due to the difficulty of gravity determination noted 
above, at least partly, or may be simply because 
our sample mostly consists of low gravity objects 
by chance. We examine if this may be due to any 
systematic effects in our UCMs. Our gravity de- 
termination is essentially based on the gravity de- 
pendence of the gas and dust opacities, but we 
are using more or less standard methods in treat- 
ing molecular and dust opacities. Also, a com- 
parative study of different modeling approaches 
did not reveal any sys tematic effect in our mod- 
els ([Helling et al.1 120081 ). Thus we cannot identify 
any systematic effect in our models. We hope to 
examine this problem further with a larger sample 
in future. 

Cushing et al. I (120081) used ultracool dwarf evo- 



lution models to transfer their (T c s,R) values to 
(T e ff,g) values, from which they obtained log 5 = 
4.7 for 2MASS J055919-1404. This result dif- 
fers from their result based on their spectral fit- 
ting referred to in Section l4~2l (log g = 5.5), but 
agrees rather well with our value given in Table [4] 
(log g = 4.5). Thus this result can be support- 
ing evidence for our log g, despite some difficul- 
ties noted above. However, the results based on 
the evolutionary sequences may suffer more or less 
similar difficulty to the spectral fitting method, 
since it depends critically on the physical parame- 
ters such as I^ff j which cannot be very accurate if 
they are based on the spectral fitting in the limited 



spectral region (as will be detailed in Section |4~51) . 

Given the radius and log g, it is in principle pos- 
sible to estimate the mass. However, SEDs are not 
so sensitive to log g and its accuracy is rather low. 
For this reason, we believe it impossible to obtain 
a reasonable mass estimates from our spectra. 

4-4-3- Critical temperature 

Inspection of Table [4] reveals that the resulting 
T cr 's show a variety of values in our limited sam- 
ple, and this result implies that the dust cloud 
properties differ significantly among different ob- 
jects. In modeling dusty dwarfs, it is still difficult 
to determine the dust cloud properties from the 
basic physics, and we had to introduce such an 
empirical parameter as T cr to represent the dust 
cloud properties, at least partly. Other models 
of dusty dwarfs also have a more or less s imilar 
feature. For example, Ackerman fc Marlevl ( 200ll ) 
introduced the cloud sedimentation efficiency / sc d 
in their models and, as discussed in Section |4.3.5[ 
their / sc( j and our T cr consistently showed that the 
dust cloud in 2MASS J055919-1404 is rather thin. 
We hope that a larger sample will be analyzed to 
clarify the nature of dust clouds and that the re- 
sult will serve as a guide towards a more physical 
theory of cloud formation in dusty dwarfs. 

4-4-4- Chemical composition 

We assumed that the chemical composition in 
such unevolved objects as brown dwarfs should be 
the same as the Solar composition. However, the 
Solar composition itself is by no means well estab- 
lished yet. For example, we adopted log Ac = 8.60 
and log A = 8.92 (on the scale of log A H = 12.0) 
in our initial version of U CM based on the latest 
results known at that time [ Tsuji 2002f ). However, 
these values were revised to be log Ac = 8.39 and 
log Ao — 8.69 based on the three-dimensional(3D) 
time-dependent hydrodynamical model of th e so- 



lar photosphere ( Allende Prieto et al. 2002; and 



our present version of UCM is bas ed on these 



revis ed values as noted previously (jTsuii et al. 
20041 ). These values are close to the more re 
cent values of log Ag = 8.43 and logAo = 8.69 



(.Wplmxi el al. 2009f ). Since a direct determina- 



tion of abundances in brown dwarfs cannot be ex- 
pected in the near future, we think it best for now 
to use the most reliable Solar composition. 
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Unlike the cases of evolved stars in which the 
surface chemical composition may suffer drastic 
variations due to convective dredge- up of the prod- 
ucts of nuclear processing in the interior, large 
variations of the surface chemical composition may 
not be expected in brown dwarfs, but a possibility 
of small variations may not be excluded. As an 
example, we examine the effect of carbon abun- 
dance on CO and C0 2 features in J083008+4828 
and J055919-1404, since CO and C0 2 features in 
these objects could not be reproduced well with 
the composition we have assumed. For this pur- 
pose, we compute two spectra with log Ac changed 
by ±0.15 dex compared to the standard composi- 
tion we have assumed, namely, log Ac = 8.39. 

The results for J083008+4828 are shown in the 
upper diagram of Figure [JJ The effect of increas- 
ing log Ac by 0.15 dex results in a considerable 
strengthening of the methane bands, but the over- 
all fitting tends to be worse (compare c with a in 
the figure). The C0 2 and CO bands at 4.2 and 
4.6 /im, respectively, however, do not show any 
strengthening. The effect of decreasing log Ac by 
0.15 dex is rather modest (see d in the figure). The 
results for J055919— 1404 are shown in the lower 
diagram of Figure [7J The effect of changing log Ac 
by ±0.15 dex can be noticed on the CH4 bands, 
but CO2 and CO bands remain almost unchanged. 
We conclude that the large discrepancies between 
the observed and predicted CO and C0 2 bands in 
J083008±4828 and J055919-1404 cannot be the 
problem of the carbon abundance. 

4-4-5- Micro-turbulent velocity 

We have assumed the micro-turbulent velocity 
to be the Solar value of 1 kms -1 throughout. Just 
to see the effect of the micro-turbulent velocity, we 
increase it to be 2 kms -1 and the resulting spectra 
are compared with those for the case of 1 kms -1 
for J083008+4828 and J055919-1404 in the up- 
per and lower diagrams, respectively, of Figure [Jj 
Clearly, dependence on the micro-turbulent veloc- 
ity is quite small (compare a and b in Figure [7J as 
expected for such high density photospheres dom- 
inated by the pressure broadening. 



4.5. Errors and limitations of the present 
analysis based on the UCMs 

Since we examine the effect of T e g by steps 
of 100 K, our T c ff cannot be more accurate than 
±50 K. Further, combined with the similar uncer- 
tainty in T cr , the final uncertainty in our estima- 
tion of T e g may be ±100 K. But we should no- 
tice that the parameters we have found are simply 
those that explain the AKARI spectra in the 2.5- 
5.0 [im region. We tried more or less similar anal- 
ysis of the 1.0-2.5 /im region of our objects and 
found some differences in the resulting parame- 
ters, although the differences are mostly within 
the errors outlined above. In other words, the so- 
lution obtained by fitting other spectral regions 
may not in general explain the 2.5-5.0 fim spectra. 
As summarized in Figure El we are trying to reach 
a unified understanding of the series of new spec- 
tra observed for the first time by AKARI covering 
from mid-L to late-T dwarfs. For this purpose, we 
believe it best to apply the parameters obtained 
from the AKARI spectra themselves and analyze 
all our sample consistently. 

As shown bv lCushing et al. ( 20081 ). the effective 
temperatures obtained by fitting a limited spec- 
tral region differ from those obtained by fitting the 
full SED (0.95-14.5 /xm) by typically « 200 K and 
by as large as 700 K in the worst case. As noted 
above, we confirm essentially the same difficulty 
with our UCM. The reason for this difficulty may 
be because the present models are far from per- 
fect. For example, we do not know the exact form 
of the dust opacities which suffer from many un- 
known effects (e.g. composition, size, shape, im- 
purity etc.) and of the molecular opacities based 
on imperfect line lists. Moreover, we do not know 
yet the exact nature of the dust clouds formed in 
the photosphere of cool dwarfs. For these reasons, 
we think it difficult to expect such rigorous numer- 
ical accuracy for brown dwarfs as realized in model 
photospheres of ordinary stars. We intended from 
the beginning that our models can at least be of 
some help as a guide in interpretation and analy- 
sis of the o bserved data of ultracool dwarfs (e.g. 
Tsuiill200lh . 

A more perfect model would allow us to more 
safely analyze a wider spectral region. However, 
current models of brown dwarf atmospheres are far 
from perfection. Fitting of a wider spectral region 
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SDSS J083008+4828 

a) Tc1700/Te1500/logg4.5/std. ucm 

b) Tc1700/Te1500/logg4.5/ Vmic = 2km/s 

c) Tc1700/Te1500/logg4.5/ logAc: + 0. 1 5dex 

d) Tc1700/Te1500/logg4.5/ logAc: -0. 1 5dex 




2MASS J055919-1404 

a) Tc1900/Te1200/logg4.5/std. ucm 

b) Tc1900/Te1200/logg4.5/ Vmic = 2km/s 

c) Tc1900/Te1200/logg4.5/ logAc: + 0. 1 5dex 
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Fig. 7. — The effects of changing log Ac by ±0. 15 dex and of increasing the micro-turbulent velocity from 1 to 
2 kms -1 are shown for J083008+4828 in the upper diagram. The same effects are shown for J055919— 1404 
in the lower diagram. 
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may sample more regions of uncertain opacities at 
the same time, and is by no means easier to have 
better results. Thus we must allow for additional 
uncertainties beyond those found by fitting one 



spectr al region. We fully agree with lCushing et al 



(|20081 ) who concluded that the accurate determi- 
nation of the physical parameters is still an elusive 
goal, even though SED fitting method shows rea- 
sonable success in general. 

For now, we must be satisfied in that our SED 
fitting analysis based on our UCM has also been 
done fairly consistently for the new spectral region 
observed by AKARI. The UCM is a semi-empirical 
model in which the property of the dust clouds is 
represented by a single parameter T cr and, for this 
very reason, our UCM is free from yet unknown 
details of the cloud formation processes. Thus, our 
UCM can be flexible enough to interpret the major 
characteristics of the observed spectra in terms of 
a few basic physical parameters. The UCM, how- 
ever, assumes LTE throughout and hence cannot 
be applied directly to the observed features show- 
ing possible deviations from LTE. We will discuss 
such a case in the next section. 
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Fig. 8. — Molecular abundances, presented as par- 
tial pressure vs total pressure (pressure in units 
of dyncm~ 2 ), in LTE (solid lines) and in non- 
equilibrium modified by vertical mixing for the 
case of b = 0.1 (dashed lines) in the photosphere 
of the T8 dwarf 2MASS J041519-0935 (Model: 
T eS /T CI / \ogg =800K/T cond /4.5). The arrow in- 
dicates the onset of convection. 
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Fig. 9. — The effect of vertical mixing on the 
CO spectrum in 2MASS J041519-0935 (T8): The 
cases of LTE, b = 0.03, and b = 0.10 are compared 
with the observed spectrum. 



5. Spectral features unexplained by the 
LTE models 

In the previous section, the overall SEDs can be 
understood reasonably well with the UCM based 
on LTE, but the observed features of CO and CO2 
in T-dwarfs (and possibly late-L dwarfs) cannot 
be fitted at all with the spectra predicted by the 
model. As for CO, this fact has been known for 



late-T dwarfs (e.g. 


Noll. Geballe. & Marlevl 1997: 


ODDenheimer et al. 


19981: Geballe et al. 12009b. 



and we confirm this fact by the better quality data 
for T8 dwarf 2MASS J041519-0935 (Figure [Bf). 
But we also find that CO bands are enhanced 
significantly already in the mid-T dwarf 2MASS 
J055919-1404 (FigureEfe) as well as in the late-L 
dwarf SDSS J083008+4828 (Figure EH). A new 
and more difficult problem is that the CO2 band 
at 4.2 /im is enhanced enormously in all the cool 
dwarfs later than L9. 

5.1. Late-T dwarf 

The unexpected detection of CO in the late- 
T dwarf Gl 229B has been interpreted as due to 
vertical mixing of CO from the warm layers where 
CO is still abundant to the cooler layers where CO 
is mostly transformed to CH4 under LTE. This is 
possible since CO is so stable that the time scale 
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Tchcm of the reaction 

CO + 3H 2 ^ CH 4 + H 2 

is ve ry slow compared to the mixing time scal e 
r mix dGriffith fc Yelle Ill999t ISaumon et al.ll2000l) . 
and a detailed computation including such a 
non-equilibrium effect has been carried out (e . g . 
Saumon et aDl2007l : iHubenv fc Burrowsll2007l ). 

We try to explain our observations by a simple 
computation as follows: We show the LTE molec- 
ular abundances of CO, H2O, CH4, and CO2 in 
the photosphere of 2MASS J041519-0935, using 
the UCM of (T eff , T cr , log 5) = (800K, T cond , 4.5), 
by the solid lines in Figure |U We assume that 
the time scale of CO destruction r c h m is equal to 
the timescale of mixing T m ; x at a depth point P 
in Figure Eland CO/CH4 abundance ratio is fixed 
at the value b of that point. We examine sev- 
eral values of b. The case of b — 0.10 is shown 
in Figure [5] as an example. The resulting non- 
equilibrium abundance of CO is indicated by the 
dashed line. If H 2 and C0 2 abundances relative 
to CO are fixed at the values at P as are CH4, 
H 2 as well as CH 4 abundances in the upper lay- 
ers should be decreased and the results are also 
shown by the dashed lines. The decreases, how- 
ever, are rather minor and difficult to see in the 
figure. On the other hand, C0 2 abundance shows 
a large increase as seen by the dashed line in Fig- 
ure M 

We compute the spectra based on the non- 
equilibrium abundances for b = 0.03 and b = 0.1 
together with that of the full LTE, and the re- 
sults are compared with the observed spectrum 
of 2MASS J041519-0935 in Figured Inspection 
of the figure reveals that the observed CO spec- 
trum of 2MASS J041519-0935 can roughly be ac- 
counted for with b = 0.1. However, the predicted 
C0 2 band is still too weak compared with the ob- 
servation even for b = 0.1. Of course, there is 
no reason why C0 2 /CH4 ratio can be fixed at the 
value at point P, since chemical reaction time-scale 
of C0 2 cannot be the same as for CO. Further, 
there is little possibility that more C0 2 can be 
dredged-up, since C0 2 abundance in the deeper 
layers is not as large as shown in Figure |5J 

On the other hand, it is possible that C0 2 
abundance will change to the equilibrium value of 
the local physical condition if the chemical time 
scale related to C0 2 formation is not as slow as 



in the case of CO. However, if C0 2 attains its 
equilibrium value with the abundant CO result- 
ing from vertical mixing, C0 2 abundance will be 
too large. In fact, C0 2 is more abundant than 
CO at logPg < 4.0 in LTE as can be inferred 
from the solid lines in Figure [8] For this rea- 
son, we cannot assume simply that C0 2 can be 
in equilibrium with CO. We cannot predict the 
precise value of C0 2 abundance without more de- 
tailed non-equilibrium analysis, but we can expect 
from the above consideration that the C0 2 abun- 
dance can be somewhat larger than those shown 
by the dashed line in Figure in the surface lay- 
ers of 2MASS J041519-0935 and hence the C0 2 
band could then be stronger than those suggested 
in Figure [9l 

5.2. Mid-T dwarf 

A new result of the AKARI observation is that 
the mid-T dwarf 2MASS J055919-1404 also shows 
much stronger CO and C0 2 bands than those pre- 
dicted with the use of the UCM based on LTE. We 
show the LTE molecular abundances in the photo- 
sphere of 2MASS J055919-1404, using the UCM 
of (T cf f, T cr , logg) = (1200 K, 1900 K, 4.5), by the 
solid lines in Figure [TOl 

We first try the vertical mixing model used for 
the late-T dwarf 2MASS J041519-0935. We again 
assume that the time scale of CO destruction T c h em 
is equal to the timescale of mixing r m ; x at a point 
P in Figure [TOl and the CO/CH4 abundance ratio 
is fixed to a value b at this point. We again test a 
few values of b. A case of b — 0.50 is shown in Fig- 
ure [TOl The resulting non-equilibrium abundance 
of CO is indicated by the dashed line. Also, H 2 
and C0 2 abundances relative to CO are fixed at 
the values at P and shown by the dashed lines in 
Figure [TO] 

We compute the spectra based on the non- 
equilibrium abundances for b = 0.1 and b — 0.5 
together with those based on the LTE values, 
and the results are compared with the observed 
spectrum of 2MASS J055919-1404 in Figure HU 
The predicted spectra for non-equilibrium CO and 
C0 2 abundances show additional depression in the 
region of CO i?+P-branches but little change in 
the 4.2 /zm region, compared to the predicted one 
for the LTE case. This result indicates that the ob- 
served CO band can be explained at least partly 
by the vertical mixing model, but the observed 
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CO2 remains completely unexplained. The ob- 
served CO feature is still too deep compared with 
the predicted one, but this may be due to the ef- 
fect of the deep CO2 feature whose origin is still 
unknown. 

We conclude that the very strong CO band ob- 
served in 2MASS J055919-1404 by AKARI can 
be at least partly explained by the vertical mixing 
model but the CO2 feature cannot at all. 
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Fig. 10. — Molecular abundances in LTE (solid 
lines) and in non-equilibrium by vertical mixing 
for the case of b — 0.5 (dashed lines) in the 
photosphere of T4.5 dwarf 2MASS J055919-1404 
(Model: (T eff , T cr , log 5) = (1200 K, 1900 K4.5)). 
The arrow indicates the onset of convection. 

5.3. Late-L dwarf 

We detected a deep absorption around 4.5 /zm 
in the late-L dwarf SDSS J083008+4828 and it 
can be identified with the CO P-branch. The P- 
branch tends to be deeper than P-branch in CO 
because of narrower line separation. We notice 
such a structure of the CO fundamental consist- 
ing of the deep P-branch and weak P-branch in 
2MASS J055919-1404, and this similarity lends 
support in identifying the 4.5 fxm feature in SDSS 
J083008+4828 as due to the P-branch of CO in 
the relatively low-quality spectrum of this source 
(see also, Figure [I}, even though the presence of 
such a deep CO feature is quite unexpected in L- 
dwarfs. We also see that the 4.2 feature is very 
strong. 
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Fig. 11. — Predicted spectra in LTE and in non- 
equilibrium by vertical mixing (b = 0.10 and 0.50) 
compared with the observed one of T4.5 dwarf 
2MASS J055919-1404. 

In SDSS J083008+4828, the observed CO band 
is quite strong while CH4 band is fairly weak. 
These observed features imply that CO is quite 
abundant, and we confirm this possibility in Fig- 
ure Q2] where LTE abundances of CO, C0 2 , CH 4 , 
and H2O are shown based on the UCM of (T s, 
T cr , logg) = (1500K, 1700K, 4.5). Inspection of 
Figure [12] reveals that the CO abundance under 
LTE is already at its maximum possible value (i.e. 
almost all the carbon is in CO). Then, vertical 
mixing, if present, cannot supply additional CO 
to the upper photosphere and has little effect on 
the CO band strengths. Thus, without any com- 
putation of the spectrum, it is clear that the verti- 
cal mixing model cannot be applied to this late-L 
dwarf. Given that an increase of the CO abun- 
dance in the photosphere can no longer be possi- 
ble (Figure [T2l , not only by vertical mixing but 
also by any other method, some changes in the 
structures of the photosphere and/or in the atmo- 
sphere are required to produce the unusually deep 
CO band observed. 

The large strengthening of the CO2 band can- 
not be due to vertical mixing either, since CO2 
abundance relative to other molecules is smaller 
in the deeper layers (Figure [T2|) and hence vertical 
mixing will reduce rather than enhance the CO2 
abundance in the upper layers. This result sug- 
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gests that a non-equilibrium process other than 
vertical mixing is needed to explain the large en- 
hancement of the CO2 band. 
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Fig. 12. — Molecular abundances in LTE (solid 
lines) in the photosphere of the L9 dwarf SDSS 
J083008+4828 (Model: (T off , T cr , log 5) = 
(1500K, 1700K, 4.5)). A sudden change of H 2 
and CO2 abundances at logP g w 6.0 is due to for- 
mation of silicate (MgSiOs). The arrow indicates 
the onset of convection. 



5.4. Vertical mixing and/or other possibil- 
ities 

Our analysis confirms that the vertical mixing 
model can be applied to T-dwarfs at least for CO. 
This result certainly shows that the vertical mix- 
ing model provides one possible way to explain 
the CO anomaly, but this fact may not necessar- 
ily imply that it is a unique solution for this phe- 
nomenon. However, our analysis cannot prove the 
applicability of the vertical mixing model to the 
late-L dwarfs. 

One problem is that the convective zone is ac- 
tually situated r ather deep i n the photosphere of 
cool dwarfs (e.g. iTsuiil 120021 ) . although it is said 
that cool dwarfs are fully convective. For exam- 
ple, the upper boundary of the convective zone is 
indicated by the arrow in Figures [8] & [TTJ1 and the 
layers where vertical mixing is expected are well 
above the convective zone. Hence vertical mix- 
ing cannot be due to convective mixing. Then, 
the problem is the unknown mechanism of efficient 



mixing in the photosphere of T-dwarfs. 

A more difficult problem is how to understand 
the enhancement of CO2 bands at 4.2 fim in late- 
ly and T dwarfs throughout. We confirm that this 
phenomenon cannot be explained by the vertical 
mixing model at all. So the question is whether 
the enhancements of CO and CO2 bands in late-L 
and T dwarfs are due to the same mechanism. It 
is still possible that an unknown process result- 
ing in the enhancement of CO2 will also result in 
the CO anomaly. In any case, the vertical mix- 
ing model cannot have general applicability, and 
we must look for a more unified model that can 
be applied to CO and CO2 in L and T-dwarfs. 
Presently, we have no solution yet as to the rea- 
son why CO2 (and also CO at least partly) are 
so strong in these objects. If non-equilibrium pro- 
cesses play major roles, we wonder why the depar- 
ture from LTE can be so significant in the brown 
dwarf photosphere, where local density is so high 
that LTE condition can be naturally satisfied. 

More recently, Stephens et al. (2009) analyzed 
the 0.8-14.5 /mi region spectra of many L-T 
dwarfs and determined not only effective temper- 
ature, surface gravity and grain sedimentation ef- 
ficiency but also verti cal gas transport efficiency , 
based on the models of lSaumon fc Marlev I (|2008f ) . 
They showed that vertical mixing improves the 
fits and the observed spectra were fitted well in 
general. It will be interesting to check whether 
their models can explain the enhanced CO bands 
in late-L dwarfs. 

6. Discussion and Concluding Remarks 

AKARI provides the first and unique opportu- 
nity to take brown dwarf spectra in the important 
wavelength range 2.5-5.0 /im continuously with- 
out interference of telluric atmosphere. In partic- 
ular, molecular bands of CO, CH4, and CO2 alto- 
gether in one spectrum present a great advantage 
for investigating chemical processes in the atmo- 
spheres of ultracool dwarfs, through quantitative 
analysis of their carbon budgets. AKARI reveals 
many remarkable features of ultracool dwarfs. Our 
observations confirm the presence of CO molecules 
in the brown dwarf atmospheres of all spectral 
types from L to very late T-type objects. The 
fact that all observations ever made for the late- 
T dwarfs detected CO molecules implies that the 
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non-equilibrium chemical composition is a com- 
mon property of the late-T dwarfs. 

It has been discussed that this non-equilibrium 
abundance is due to vertical mixing in the upper 
photosphere. The CO fundamental band in the T8 
dwarf 2MASS J041519-0935 and the T4.5 dwarf 
2MASS J055191-1404 can be accounted for by 
the vertical mixing m odel (Section | 5. I I). The idea 



is also supported by Saumon et al. ~ i|2006l ) . and 



by IMainzer et al. ( 20071 ) based on the analysis of 



the mid-infrared spectrum of T7.5 dwarf Gl 570D 
taken by Spitzer /IRS. They found that the abun- 
dance of the NH3 molecule in the object is almost 
one order of magnitude smaller than that expected 
under equilibrium conditions. It is considered that 
fresh N2 molecules are continuously provided from 
inside. Also, the vertical mixing model is shown to 
be consistent with the photometry of a larger sam- 
ple of brown dwarfs obtained with Spit zer /IRAC 
(Pat ten et al. I l2006l ILeggett et al.|[2007h . 

However, our observations reveal that the CO 
fundamental band is stronger than predicted 
by the LTE models in the late-L dwarf (SDSS 
J083008+4828). Since, unlike in T-dwarfs, CO 
is already highly abundant in the upper photo- 
sphere of these objects, no appreciable increase of 
CO abundance and hence of the CO band strength 
can be expected by additional CO from the deeper 
warm layers. For this reason, some mechanism(s) 
other than vertical mixing should be looked for in 
late-L dwarfs. This quest is further strengthened 
by the detection of CO2 by AKARI as discussed 
below. 

A new piece of information provided by AKARI 
is the detection of CO2 molecule. The CO2 band is 
never observable from the ground, and its 4.2 ^m 
band is just outside of Spitzer /IRS coverage. Gas 
phase CO2 in stellar atmospheres was first de- 
tected in late-type giants by ISO /SWS observa- 
tions. It was rather unexpected that in addi- 
tion to the stretching-mode absorption band at 
4.2 /jja, the bending mode ro-vibrational bands 
in the 15 /im region are oft en seen in these 
stars, sometimes in emission ( Ryde et al. 19981: 



ljusttanont et al. Il998h . The CO2 band is gen- 
erally stronger than those expected from the 
LTE models, implying that some kinds of non- 
equilibrium processes, either formation of the 
molecule or radiative excitation, might play a 
role in the atmospheres of such stars. The sit- 



uation is of course quite different in brown dwarfs, 
with lower temperature and much higher surface 
density than in red giants. Some amount of the 
CO2 molecule is expected under LTE conditions 
in brown dwarf atmospheres. Though the relative 
abundance with respect to that of H2O and CH4 is 
small (FiguresElCroi &[JJ), the high /-value of this 
molecule helps to form a rather strong absorption 
band at 4.2 /im region as shown in Section |4~T1 (see 
Figure [5]). However, a simple comparison with the 
model tells us that observed absorption is much 
stronger than expected. 

We estimate the temperatures and column 
densities of the extra CO2 and CO absorption 
components in 2MASS J055919-1404 and SDSS 
J083008+4828. The observed spectra are divided 
by the best fit UCM spectra shown in Figure [6j 
and the residual spectra are fitted by two-layers 
plane-parallel model with given excitation temper- 
ature T cx and column density N for each molecule. 
In this simple model we adopt a Gaussian line 
width of FWHM = 30 kms -1 based upon mea- 
surements of H2O lines by high-resolution spec- 
troscopy in the J-band ([McLean et al.l 120071 ). Be- 
cause of the simplicity of the model, the fitting 
results are suitable mainly for order of magni- 
tude discussions. The excitation temperature 
is only mildly constrained by the band shape. 
Nevertheless, we obtain (T ex , N)qo 2 = (800 
K, 5 x 10 17 cm" 2 ) and (T ox , N) co = (1000 K, 
2 x 10 20 cm- 2 ) for 2MASS J055919-1404, and 
(T ex , iV) C o 2 = (800 K, 2.5 x 10 17 cm- 2 ) and 
(T ex , iV) co = (1000 K, 1.5 x 10 19 cm- 2 ) for 
SDSS J083008+4828, respectively. The results 
for 2MASS J055919-1404 are presented in Fig- 
ure H3] as an example. Considering that CO 
molecules with higher excitation temperature 
may located in a slightly interior portion of the 
photosphere, the ratio of the column densities 
N C o 2 / N co = 2 x 10~ 2 - 3 x 10~ 3 indicates lower 
limits of the abundance ratio of the two molecules. 
It is found that CO2 is obviously overabundant 
compared to CO than predicted by the UCM (Fig- 
ures [10] & [12]) . This implies that non-equilibrium 
processes other than vertical mixing may play a 
role in determining the molecular abundances es- 
pecially for C02- 

Another possibility is the presence of extended 
atmosphere beyond the hydrostatic photospheres, 
similarly in the case of red-giants. It is interest- 
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ing that Spitzer /IRS observation did not detect 
any clear CO2 bands in the 15 /im region in any 



brown dwarfs ( Rocllig et al.l l2004t iC ushing et al.. 
2006; ISaumon et alj|2006t iMainzer et al.ll2007t) . If 



the CO2 molecular layer is extended, emission be- 
yond the photosphere could compensate for the 
absorption and the observed band strength would 
be weakened. This effect is usually more effective 
at longer wavelengths. Of course, there is a more 
conventional explanation that the heavy opacity of 
H2O veils the wavelength region completely. The 
H2O opacity reaches a minimum around 4 /im and 
makes the 4.2 /im band visible. If this is the case, 
the CO2 molecules would be located in the same 
layer, or in layers internal to the H2O. This may 
be a constraint on constructing chemical models 
of the brown dwarf atmospheres. High dispersion 
spectroscopy in the 4 & 15 /im region from space 
is an ideal tool for investigating this problem. 
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Fig. 13.— The AKARI spectrum of 2MASS 
J055919-1404 is divided by the best fit UCM 
model and shown in black. Blue line indicates 
the fitted model spectrum by a plane-parallel con- 
figuration with the excitation temperature and 
the column density of (T G x, N)co 2 = (800 K, 
5 x 10 17 cm- 2 ) and (T cx , iV) co = (1000 K, 2 x 
10 20 cm -2 ), respectively. 

Although non-equilibrium conditions for the 



CO/CH4 and NH3/N2 ratios appeared to be ex- 
plained by the vertical mixing model at least in 
late-T dwarfs, the unexpected detection of CO2 
by AKARI casts doubt as to whether the verti- 
cal mixing model alone could explain the carbon 
chemistry even in late-T dwarfs. In fact, there is 
no reason why we expect excess CO2 by the ver- 
tical mixing model (Section [52]) and there should 
be some other mechanism to produce excess C02- 
A yet unknown mechanism may affect not only 
CO2, but also CO as well. In fact, such a new 
mechanism could also be what we need to explain 
the enhancement of CO in late-L to mid-T dwarfs 
for which the vertical mixing model failed. 

We continue the AKARI brown dwarf pro- 
gramme with the near-infrared channel of the IRC 
in the "post-Helium" phase (Phase 3). More 
targets are being observed to cover the spectral 
ranges from M, L, to T. With the enhanced data 
set we hope to better understand the some of the 
current mysteries of ultracool dwarfs. 
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